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Abstract 


A  theoretical  and  computational  analysis  of  the  quantum  dynamics  of  charge  carriers  in 
presence  of  electron-phonon  interaction  based  on  the  Wigner  function  is.  here  applied  to  the 
study  of  transport  in  mesoscopic  systems.  Numerical  applications  are  shown  for  a)  a  wave 
packet  .scattering  with  phonons  while  crossing  a  potential  profile  and  b)  electrons  scattering 
with  phonons  in  a  finite  device  with  open  boundary  conditions. 
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■  ditions 


2 


1  Introduction 

A  rigorous  quantum  transport  theory  of  electrons  in  semiconductor  structures  including  both 
coherent  propagation  and  scattering  mechanisms  is  still  lacking  even  though  its  need  has  been 
recognized  since  many  years.  In  particular  we  miss  a  theory  feasible  for  numerical  calculations 
to  be  compared  with  experimental  results  and/or  with  results  of  the  semiclassical  theory  with 
the  purpose  of  enlightening  the  different  predictions  of  the  two  models. 

In  the  literature  different  approches  have  been  used  [1,  2,  3].  The  Wigner  function  (WF)  T] 
provides  a  rigorous  quantum-mechanical  tool  and  constitues  a  direct  link  between  quantum 
and  classical  descriptions  of  the  evolution  of  the  system  in  phase.  In  this  approach  however 
electron-phonon  interaction  has  been  treated  using  either  relaxation  time  [4,  5]  or  a  classical 
Boltzmann  collision  operator  [6].  In  this  paper  several  advancements  in  the  efford  to  solve  the 
above  problem  exactly  are  presented. 

The  definition  of  the  WF  can  be  extended  to  include  electrons  interacting  with  phonons  as 
follows  [7]: 

fwir,p,ng,n'^,t)  = j  dr'e~^^’''p{r  +  r/2,ng;r-r'/2,n'^)  (1) 

where  r  and  p  are  the  electron  position  and  momentum,  p  is  the  density  matrix  for  the  system 
formed  by  an  electron  and  the  phonon  gas,  and  riq  is  the  occupation  number  of  the  phonon  mode ' 
q.  Trace  over  the  phouon  variables  will  lead  to  the  traditional  electron  WF. 

In  this  paper  we  discuss  an  integral  equation  for  the  above  WF  and  apply  it  to  study  quantum 
electron  transport  in  open  systems  in  presence  of  scattering.  In  particular  suitable  boundary , 
conditions  can  be  used  which  restrict  the  analysis  of  the  WF  inside  a  finite  domain.  Wigner 
trajectories  have  been  identified  also  for  the  case  of  electron-phonon  interaction  and  constitute 
a  guiding  criterion  for  replacing  initial  conditions  with  boundary  conditions  for  the  integral 
equation. 
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2  Integral  Equation 


Let  us  consider  a  system  described  by  the  Hamiltonian 

H  =  Ho+H'  (2) 


where  Ho  is  the  “unperturbed”  Hamiltonian,  containing  the  description  of  the  free  electrons, 
free  phonons  and,  in  the  present  approach,  also  the  potential  profile  inside  the  system  (device). 
H'  contains  the  scattering  agents,  i.e.  the  electron-phonon  interaction. 

The  time  evolution  of  the  density  matrix  is  described  by  the  Liouville-von  Neumann  equation. 
Starting  from  this  equation  we  derived  the  following  integral  equation  for  the  WF: 


Ur,  p,  n„  <,  t)  =  fnn’ir,  f  dr'  [  dp'Uiv',  p') 

'xfw{r',p',ng,n'g,to)  +  h^J2frin'ir,p)  ( 

^  XI  [  ^^'  [  <ip' {'^'inng,rnTng)fUir',p')Ur',.p',rnq,n'g,t') 

mrriq 

fnrrSj-  >  P  5  P  1  ^<7?  ^(7?  ^  ^  (3) 


In  the  above  equation  K'  —  and  are  its  matrix  elements  on  the  basis 

>}  of  Ho.  The  coefficients  fimi^^p)  are  given  by 


If  p  /  r^ 


> 


(4) 


They  allow  to  move  from  the  WF  to  the  density- matrix  representation  and  viceversa.  Finally 
/ia;(n,ng).  indicates  the  total  unperturbed  energy  of  the  electron-phonon  system  when  the  elec¬ 
tron  is  in  the  n-th  eigenstate  and  the  phonon  bath  is  in  a  state  with  n  phonons  in  mode  q. 

In  Eq.  3  the  first  term  in  the  r.h.s.  describes  the  ballistic  coherent  propagation  from  the 
initial  time  to  to  the  observation  time  t.  The  second  term  describes  the  contribution  of  the 
unknown  WF  to  the  ‘Hast”  interaction  vertex  at  any  time  t'  between  to  and  t  followed  by  the 
free  propagation  from  t'  to  t.  Even  though  the  equation  is  written  for  the  entire  system  and 
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contains  also  the  phonon  variables,  it.  has  been  shown  that  a  Monte  Carlo  solution  allows  for  a 
proper  average  over  the  phonon  coordinates,  as  long  as  hot-phonon  effects  arc  neglected  [8]. 


3  Wigner  Trajectories 


It  is  well  known  that  a  point-like  electronic  WF  in  phase  space  cannot  represent  a  real  physical 
system,  since  it  violates  the  uncertainty  principle.  We  may  however  consider  the  evolution  in 
time  of  a  (5— like  contribution  to  the  WF  inside  Eq.  3.  For  free  electrons  its  trajectory  in  phase 
space  is  the  same  as  for  semiclassical  particles.  In  fact,  for  a  plane- wave  basis  the  coefficients 
fim  are  given  by 

fkk'{r,p)  =  +  k'))  (5) 

and  the  ballistic  evolution  of  the  WF  as  given  by  the  first  term  on  the  rhs  of  Eq.3  yields: 


fw{r,p,ng,ng,t)  =  fyj{r- ^{t-to),p,ng,n'^,to)e  “K)){* 


(6) 


where  m  is  the  electron  effective  mass.  The  first  factor  and  the  exponential  factor  describe  the 
free  trajectory  of  the  electron  and  the  time  evolution  of  the  free-phonon  bath,  respectively. 


This  result  is  not  surprising  since  in  absence  of  scattering  and  up  to  quadratic  potentials  (har¬ 
monic  oscillator)  a  differential  equation  can  be  written  for  the  WF  that  coincides  with  the 
Boltzmann  equation,  so  that  the  dynamical  evolution  of  the  WF  follows  the  same  trajectories  of 
the  representative  points  of  a  classical  gas.  This  implies  a  deformation  of  the  WF  while  evolving 
in  time  due  to  the  fact  that,  as  time  increases,  the  higher-momentum  components  move  faster 
than  the  lower- momentum  components. 


If  phonon  scattering  is  included,  similar  trajectories  exist,  with  very  interesting  properties,  as 
long  as  a  single  interaction  diagram  is  considered  without  time  integrations.  In  order  to  prove 
this  result  Eq.  3  must  be  substituted  into  itself  to  obtain  the  second-order  term  corresponding  to 
one  scattering  event.  Then  assuming  a  constant  potential  and  using  again  plane  waves  as  basis 
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set,  the  second-order  correction,  e.g.  for  the  terms  corresponding  to  phonon  emission,  results  to 
be  - 

=  2^lj2.F\q)  <  ji,  >  f' dt'  f  df 

V  g  ^  to  J  to 

here  for  simpLcity  a,  one-dimensional  case  is  considered  and  F{q)  is  the  coupling  function 
entering  the  electron-phonon  hamiltonian  [9].  The  trajectories  described  by  the  above  equation 
are  represented  in  the  diagrams  contained  in  Fig.  1.  They  correspond  to  the  semiclassical 
trajectories  where  g/2  is  transferred  to  the  electron  at  each  vertex.  As  a  result  of  the  time 
integrations,  hov.ever,  a  sum  over  infinite  trajectories  gives  the  WF  correction  corresponding  to 
a  single  scattering  event. 

4  Initial  and  Boundary  Conditions 


Guided  by  the  considerations  of  the  previous  section  we  face  the  problem  of  finding  an  equation 
equivalent  to  Eq.  3  for  the  WF  in  a  finite  region  inside  given  boundaries.  The  problem  is 
not  trivial  owing  to  the  nonlocality  of  the  interaction  and  to  the  integration  from  -oo  to  -fee 
in  the  integral  over  z'.  We  have  proved  that  for  the  solution  of  the  integral -equation  it  is 
possible  to  substitute  the  knowledge  of  the  WF  over  all  space  at  a  given  (initial)  time  to  with 
the  knowledge  of  the  same  function  inside  the  region  of  interest  at  to  and  at  all  times  at  the 
boundary  for  “entering”  momenta.  The  integral  equation  for  such  a  situation  results  to  be 


/.u  = 


+  /l 


Ylir-.'.'{z,Pz)  [  dt'< 

Tin* 


■~i[uj{nnq)-uj(^n' n'^)){t-t' ) 


r+oo 


cA-  -  f  dp',  (A,p'J  A  (A,p'„nq,n'  t') 

'  J  — OO  ^  ^ 
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+  /  /  ^p'z  ,  rnmg )  {z' ,  p'^ )  /„  (2' ,  p, ,  m, ,  n', ,  t')  ' 

f nm  ('^  ^  P2;)  fw  5  '^Q'i  t  ^  1~L  ^TTlTTlq^  Tl  ]  ^  ' 

where  again  for  simplicity  a  one-dimensional  case  has  been  considered,  and  —A  and  A  are  the 
device  boundaries. 

5  Applications 

As  a  case  study  we  have  first  analysed  a  single  electron-phonon  scattering  event  for  a  wave 
packet  crossing  a  potential  profile  [10],  In  this  case  the  numerical  procedure  is  simplified  by 
the  fact  that  the  whole  system  can  be  enclosed  in  a  large,  but  finite  box.  Thus  the  basis  set 
of  functions  is  discrete  and  can  be  used  without  finite-difference  approximations  in  k-space. 
Three  cases  have  been  considered  for  a  single  electron-phonon  scattering  event:  a)  the  effect  of 
a  quantum  collision  over  the  free  propagation,  b)  intra-collisional  field  effect  when  the  phonon 
interaction  takes  place  during  a  propagation  through  a  region  with  a  constant  electric  field, 
and  c)  scattering  during  the  otherwise  resonant  propagation  through  a  double  barrier.  As  an 
example  Fig.  2  shows  the  scattering  and  scattered  WF  at  t  =  l.bps  after  the  initial  condition 
for  an  electron  crossing  a  double  barrier  in  resonant  conditions  and  emitting  an  optical  phonon. 
The  interaction  is  switched  on  when  the  packet  hits  the  double  barrier,  namely  at  t  —  0.5ps 
after  the  initial  condition.  A  complex  quantum-dynamical  interference  between  potential  profile 
and  scattering  is  observed  at  successive  times  [10]. 

In  order  to  verify  the  correctness  of  the  use  of  the  boundary  conditions  for  an  open  mesoscopic 
system  as  described  in  Section  4  we  have  first  evaluated  the  ballistic  evolution  of  the  first  term 
in  Eq.8  assuming  a  step  potential  profile,  /^  =  0  at  t  =  to  inside  the  device,  and  a  Maxwellian 
local-equilibrium  incoming  distribution  at  the  left  and  right  boundaries.  Results  are  showm  in 
Fig.  3  at  different  times  after  the  initial  condition. 
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At  the  longest  times  (not  shown  in  the  figure  for  space  reasons)  we  recover  for  the  same  result 
which  is  obtained  without  using  boundary  conditions  for  a  stationary  infinite  system,  where 
is  evaluated  by  means  of  an  equilibrium  density  matrix  diagonal  over  the  scattering  states. 

In  order  to  describe  quantum  transport  across  a  very  small  device  in  presence  of  scattering 
the  integral  Eq.8  is  expanded  perturbativejy  through  iterative  substitutions.  The  second-order 
contribution  is  the  sum  of  four  terms:  real  and  virtual  absorptions  and  emissions.  As  an  example 
Fig.  4  shows  the  result  for  the  in-emission  term  summed  over  all  phonon  modes  q  which  compares 
very  well  with  the  classical  result. 

An  interesting  difference  instead  emerges  when  only  a  single  mode  q  is  considered  in  the  quantum 
calculation,  as  shown  in  Fig.  5.  At  short  times  right  after  the  entrance  energy  non-conservation 
modifies  the  in-scattering  term. 


Partially  supported  by  ARO  and  ONR  through  ERO  through  Contract  n.  N68171-96-C-9089. 
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Figure  1:  Example  of  the  VVigner  trajectories  that  contribute  to  the  result  in  Eq.  3  for  real  and 
virtual  emissions  of  a  mode  q. 


Figure  2:  Unperturbed  Wigner  function  (above)  and  second-order  correction  due  to  a  quantum 
phonon  emission  (below)  for  an  electron  crossing  a  double-barrier  potential  profile  in  resonant 
conditions  at  different  times.  The  interaction  is  active  from  t  =  0.5  ps  after,  the  initial  condition. 


Figure  3:  Ballistic  evolution  of  the  first  term  in  Eq.  (8)  at  d  times  (t  =  10  fs  and  8O/5)  for 
electrons  entering  from  the  boundaries  into  a  region  with  a  potential  step. 


Figure  4:  In-emission  term  summed  over  all  phonon  modes  g  at  ^  ~  50  fs  after  the  electron- 
phonon  interaction  has  been  switched  on. 


Figure  5:  Contribution  to  the  Wigner  function  due  to  real  emission  of  a  mode  q  =  1.66  •  10^ 
at  200/s  after  the  electron-phonon  interaction  has  been  switched  on.  Only  electrons  entering 
with  positive  p  are  able  to  emit.  At  long  times  classical  conservation  is  recovered. 
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